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Using low temperature scanning tunneling spectroscopy, we probe the Landau levels of the topologically 
protected state of Sb2Te3(0001) after in-situ cleavage of a single crystal. Landau levels are visible for magnetic 
fields B > 2 T at energies, which confirm the Dirac type dispersion including the zeroth Landau level. We 
find different Dirac velocities for the lower and the upper part of the Dirac cone in reasonable agreement with 
previous density functional theory data. The Dirac point deduced from the zeroth Landau level shifts by about 
40 meV between different areas of the sample indicating long range potential fluctuations. The local potentials 
are correlated to different local defect densities varying slightly stronger than expected from a statistical distri¬ 
bution. Moreover, the width of the Landau level peaks is analyzed. It is found to increase, mostly linearly, with 
the energy distance to the Fermi level. Consequently, we attribute the peak width to a dominating scattering of 


the hot quasiparticles by electron-electron interaction. 

PACS numbers: 71.20.Nr, 71.70.Ej, 73.20.At 

I. INTRODUCTION 

Topological indices are a new paradigm to classify solids 
by relating bulk properties unequivocally to the conductivity 
at the rim of the sample They have been used to cat¬ 

egorize the quantum Hall effect 0111 as well as to predict a 
non-magnetic quantized transversal conductance in superflu¬ 
ids 0. More recently, they led to the experimental discovery 
of two-dimensional (2D) topological insulators (TIs) EHm, 
strong and weak three-dimensional (3D) TIs lfT2UT7ll . topo¬ 
logical crystalline insulators mini and the anomalous quan¬ 
tum Hall effect EOl 1211 . The strong 3DTIs are so far the 
most versatile class in terms of different realizations in ma¬ 
terials ll3H5l l^l^ . The driving force within these materials 
is a strong spin-orbit (SO) interaction leading to a partial in¬ 
version of bands around the band gap, while respecting the 
time-reversal symmetry. As a consequence, non-trivial sur¬ 
face states emerge within the bulk energy gap, forming an odd 
number of Dirac cones and exhibiting a chiral relationship be¬ 
tween spin and momentum lfT2l l24l - E7l . The Landau quan¬ 
tization of the massless Dirac electrons in a magnetic field 
B exhibits a square-root dependence with respect to B and 
a held-independent zeroth Landau level (LL) ll28l l29ll . both, 
in contrast to gapped two-dimensional electron systems with 
parabolic dispersion. In that sense, the zeroth Landau level 
called LLo is a hngerprint of the topological protection of the 
surface states. 

Low-temperature scanning tunneling spectroscopy (STS) can 
map the Landau quantization down to the atomic scale ex¬ 
hibiting discrete energy peaks in the differential conductiv¬ 
ity (dl/dV) spectra. This method has been applied, hrstly, 
to semiconductor systems with parabolic dispersion Il30l43^ . 
and later to Dirac electron systems, such as graphene imiiTi- 
1^ . strong 3DTIs 3D topological crystalline insu¬ 

lators ll45l l46l and 3D Dirac semimetals BtII . The method 
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has been used to determine the energy dispersion including its 
lifetime broadening |l28l[30l|3ll|32]||8]|43|4T]|43l|45}i7l, to 
probe potential fluctuations ll^[Mll42ll . to probe the Landau 
level wave functions MM, to study the related quantum 
Hall effect on the local scale 13^ O, and to observe influ¬ 
ences of electron-electron interaction ll3^ . 

In this work, we apply the method to the phase change 
alloy Sb 2 Te 3 BSll . being a prototype 3DTI with only one 
spin-polarized Dirac cone located at the L-point, as has been 
revealed by spin- and angle-resolved photoemission spec¬ 
troscopy (spin-ARPES) |l49] |50l, conventional ARPES lISTl 
l52l and the de-Haas-van-Alphen effect ll53l . Similar to other 
binary TIs, Sb 2 Te 3 has an intrinsic defect population which 
pins the Eermi level Ep within the bulk valence band. Com¬ 
bined STM and density functional theory (DFT) calcula¬ 
tions of Sb 2 Te 3 have revealed Sb vacancies (Vsb) and Sb-on- 
Te antisites (SbTo) as the energetically favorable p-type de¬ 
fects. These defects appear, e.g., after molecular beam epitaxy 
(MBE) of Sb 2 Te 3 thin hlms iMl and are responsible for the 
natural p-type conductivity 

Previous STS data on MBE grown Sb 2 Te 3 thin films of 
7 quintuple layers (QLs) Il43l have found LLs corroborat¬ 
ing the Dirac cone nature of the topological surface state, 
where the overall Dirac velocity was determined to be ud = 
4.3 • 10® m/s. Moreover, the energy dependence of the LL 
width was found to be compatible with a dominating electron- 
electron scattering ||43| . Here, we show, that despite a higher 
defect density in bulk Sb 2 Te 3 of (2 — 4) • 10®^ cm“^, STS 
reveals Landau quantization already at a magnetic held of 
B > 2 T in bulk samples. The LLs show a similar Dirac- 
like dispersion of the topological surface state (TSS) as found 
in l43l . Additionally, we reveal different Dirac velocities tiD 
for the lower and the upper part of the Dirac cone. The distinct 
riD values differ by about 20 %, and are both in good agree¬ 
ment with previous DFT calculations and ARPES data of the 
hole part of the Dirac cone of the same bulk crystal B9l . The 
Dirac point energy Ed deduced from the energy of LLq spa¬ 
tially varies by up to 40meV with respect to Ep, which could 
be directly related to different amounts of local p- and n-type 
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FIG. 1. (color online) (a) Large scale STM image of cleaved 
Sb2Te3(0001) (V — 0.9 V, I — 50 pA) revealing large terraces 
of widths up to 50 nm. Inset: Line profile along the blue line of 
the main image showing the step heights which correspond to the 
height of one QL (~ 1 nm). (b) Atomically resolved STM image 
(V = 0.4 V, I = 1 nA) recorded in the area marked by the black box 
in (a). A hexagonally arranged pattern of the top surface Te atoms is 
observed with an average atomic distance of 0.42 nm. Different types 
of defects are visible as clover-shaped structures appearing dark and 
bright. 


defect densities. The quasiparticle lifetime as deduced from 
the peak widths AE of the LLs reveals an approximately lin¬ 
ear dependence of AE on energy with respect to Ep. This 
is attributed to a dominant contribution of electron-electron 
interaction to the inelastic scattering rate by a detailed anal¬ 
ysis of the expected features of different scattering channels, 
thus, generalizing the conclusion of ||43]| towards larger defect 
densities and, hence, putting it on more solid grounds. 


II. EXPERIMENT 

Scanning tunneling microscopy (STM) and STS measure¬ 
ments are performed in a home-built microscope in ultra-high 
vacuum (UHV) at a temperature of T = 6 K ll56ll and a variable 
B field perpendicular to the sample surface up to i? = 7T. 
The Sb 2 Te 3 single crystals are cleaved in UHV, at a base 
pressure of 1 • 10“^°mbar, and are transferred into the pre¬ 
cooled STM directly afterwards. STM topography images 
are recorded in constant-current mode at current / and sample 
voltage V. The local density of states (LDOS) of the sample 
surface is measured via STS by locally resolved dI/dV{V) 
curves using lock-in technique with modulation frequency 
/mod =1.5 kHz and amplitude Vnod = 2-4 mV„is- This results 
in an energy resolution of JE « (3, ikp,TY + (1, SeUmod)^ 

« 4-7 meV E?) as has been crosschecked for the used STM 
previously ll5^ . The dl/dV{V) spectra are obtained by stabi¬ 
lizing the tungsten tip at a distinct tip-sample distance defined 
by the current /stab and sample voltage Vtab- 


III. IDENTIEICATION OE DEEECTS 


A large-scale STM image of the cleaved Sb2Te3(0001) sur¬ 
face is shown in Fig. [TJa). It reveals typical terraces with 
widths of up to 50 nm separated by steps of « 1 nm (inset of 
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FIG. 2. (color online) (a) dl/dV{V) spectrum recorded on cleaved 
Sb2Te3(0001) at B = 6.7 T showing Landau quantization of the 
topological surface state (Utah = 0.3 V, /stab = 400 pA, LVod = 
4 mV). The spectrum is an average of ten consecutive spectra mea¬ 
sured at the same point on an area away from defects. Eu marks 
the position of the Dirac point located at LLq, the numbers mark the 
LL index n, and the dashed lines mark the peak positions deduced 
from Lorentzian fits. Inset: dI/dV{V) spectrum at 73 = 0 T mea¬ 
sured at the same position as the one in the main image with Eo 
marking the minimum of the curve (Utah = 0.3 V, /stab = 50 pA, 
Tmod = 4 mV), (b) Landau level energies deduced from Lorentzian 
fits of the peaks in (a) and plotted against sgn{n)^/\n\B. The line is 
a linear fit resulting in the Dirac velocity vd as indicated. 
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Fig. [2 a)), which corresponds to the height of one QL. More¬ 
over, the atomic structure of the surface is visible at smaller 
scale (Fig. [TJb)) and reveals a hexagonally arranged pattern 
of Te-atoms with an atomic distance of a = 0.42 nm. Typi¬ 
cal intrinsic defects of the top QL are visible as clover-shaped 
darker and brighter areas in the STM image. The bright de¬ 
fect corresponds to a Sb-on-Te antisite (SbTe) in the Te sur¬ 
face layer, while the dark defect is a Sb vacancy located in the 
underlying Sb-layer (Vsbi). as has both been previously iden¬ 
tified by Jiang et al. Il54l . These particular defects have been 
found to be responsible for the natural p-type conductivity of 

Sb2Te3 lai. 
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IV. DEFECT DENSITIES COMPARED WITH DIRAC 
POINT ENERGIES AS DEDUCED FROM LL 
SPECTROSCOPY 

The electronic structure of Sb2Te3(0001) is, firstly, probed 
at B = 0 T via the dI/dV{V) spectrum shown in the in¬ 
set of Fig.j^a). The minimum in the differential conductivity 
at 170 meV above E-p is attributed to the Dirac point energy 
Ep, through comparison with the energy dispersion of the TSS 
from DFT calculation ||49|. In the same way, the finite inten¬ 
sity for energies below the Dirac point can be attributed to the 
bulk valence band. Thus, the STS spectrum reveals a signif¬ 
icant hole doping in accordance with other STM and ARPES 
studies Il43l[50ti52ll5^l55]l . The related large carrier density 
of 2 - 5 • 10i®cm-3 ISSlESl in combination with the sur¬ 
face Dirac cone favorably suppresses tip induced band bend¬ 
ing ll5^ . However, the found Ep, — Ep deviates from a previ¬ 
ous ARPES study on the same Sb 2 Te 3 crystal for an unknown 
reason ll49ll . 

The Dirac fermion nature of the TSS can be probed by LL 
spectroscopy, exploiting that the electron energy En in a per¬ 
pendicular i3-field is quantized into discrete values according 

to 1321 SOI ST] IMl 

En = Ed + sgn{n)hvp^2e\n\B/h, n = 0, ±1, ±2,... 

( 1 ) 

with n being the Landau level index. This equation includes 
the 73—field independent energy of LLg at n = 0, which is 
protected by the Berry phase of tt of the TSS. Eigure |^a) 
shows the dl/dV{V) spectrum measured at 73 = 6.7 T and at 
an area away from defects. Note, that the LL spectra measured 
on top of a defect closeby are barely changed. The quantiza¬ 
tion into LLs appears as peaks with non-equal energy spacing. 
The latter reveals the non-parabolic dispersion of the Dirac 
cone. The central peak is found at about V = 166 mV, which 
is exactly the energy position of Ep deduced from the dl/dV 
curve at 73 = 0 T (inset of Eig. |^a)). This peak is, thus, 
identified as LLg. Below LLg, only three further peaks are 
resolved probably due to the overlap of the lower part of the 
Dirac cone with the bulk valence band ||49|. This is in quali¬ 
tative agreement with a recent ARPES work, which proposes 
that a change in the TSS wave function occurs at the crossing 
point between the TSS and the bulk valence band due to hy¬ 
bridization M- The character of the TSS wave function is 
then expected to become more bulk-like leading to less pro¬ 
nounced LLs ll62ll . Above LLg, up to six peaks, not present at 
73 = 0 T, are discernible. 

The dispersion of the probed electron states can be de¬ 
duced by plotting the found LL energies E^ with respect to 
sgn(n)v/|n|B (eq. 0- Therefore, the dl/dV curve of Fig. 
I^a) is fitted by a Lorentzian function for each peak leav¬ 
ing the peak positions and the peak widths as free parame¬ 
ters. The resulting peak positions 73„ are indicated by dashed 
lines in Fig. j^a). The respective plot 73„(sgn(n)-\/|«|73) 
is shown in Fig. [^b). For the Landau levels around LLg, it 
reveals the expected linear dispersion close to Ep confirm¬ 
ing the massless Dirac fermion character of the TSS. Towards 
higher LL index (n > 4, n < -3), the dispersion gets slightly 



FIG. 3. (color online) (a) STM image (V = 0.9 V, I = 100pA) re¬ 
vealing typical types of defects at the surface of Sb 2 Te 3 ( 0001 ). (b) 
Zoom into the marked area of (a) showing a close up view with four 
different types of defects marked by different colored ellipses, which 
are labeled according to IMI 


non-linear. A non-linearity away from T^d has also been found 
in previous DFT calculations of Sb2Te3(0001) ll49l . The re¬ 
sulting vp deduced from a linear fit to the data amounts to 
ud = 4.2 ± 0.1 • 10^ m/s, which nicely agrees with the value 
obtained by LL spectroscopy on a 7 QL film of Sb 2 Te 3 grown 
by MBE (UD = 4.3 • m/s) iH. 

Next, we determine the local defect density within the top 
QL. This defect density will locally determine the electrostatic 
potential as measurable by 73 d — Ep with Ep being the en¬ 
ergy of LLg. We rely on the comparative STM and DFT study 
mentioned above ll43l . This study has identified the Sb va¬ 
cancies (Vsb«^) in the two different Sb-layers within the top 
QL (n = 1,2: layer label), and the Srixc antisites located at 
the top Te-layer. These defects are acceptors, while the also 
identified Te-on-Sb antisites (Tesb) act as donors. Note, that 
Sbxe antisites on subsurface Te layers of the top QL are not 
identified within our STM data, even though they should be 
visible with reduced intensity at the surface, most probably 
due to their higher formation energy ll43l . 
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TABLE I. Number of specific p- and n-type defects (as marked in Fig. [^b)) found in six distinct areas, all being 40 x 40 nm^ in size. The 
counting uncertainty amounts to 10 %. The overall defect density tid is the sum of all defects divided hy the prohed area. The effective p-type 
doping Peff is the sum of the p-type defects minus the number of the n-type defects divided by the product of probed area and thickness of 
a QL (1 nm). The energy position of the Dirac point Eu for area 1 and 2 is deduced from the fitted position of LLq. Area 1 (Area 2) is the 
region of the LLs shown in Fig.|^a) (Fig.|^a)). 



Area 1 

(LLs in Fig.[2|a)) 

Area 2 

(LLs inFig.|4l;a)) 

Area 3 

(area within Fig.|3|a)) 

Area 4 

Area 5 

Area 6 

, ^ Vsh 

32 

56 

48 

53 

49 

32 

p-type detects: 

oUTe 

14 

6 

10 

3 

3 

8 

n-type defects: Tesb 

3 

3 

6 

2 

4 

2 

overall defect density 

0.031 nm^^ 

0.041 nm”"^ 

0.040 nm~^ 

0.036 nm^^ 

0.035 nm“^ 

0.026 nm^^ 

effective p-type doping 

2.7 • 10"^ cm^"* 

3.7 • 10"'^ cm^^ 

3.2 • lO^^'cm’^ 

3.4 • lO^^'cm^^ 

3.0 • 10"^ cm^^ 

2.4 • lO^'^cm”^ 

1 

a 

(166.1 ± 0.05) meV 

(204 ± 1) meV 

/ 

/ 

/ 

/ 


Figure[^a) shows a STM image of the Sb2Te3(0001) surface, 
resolving different types and numbers of clover-shaped de¬ 
fects. The apparent defect density is 3.9 • 10^^ cm“^, which is 
slightly larger than defect densities found in MBE grown thin 
films of SbsTea lEl (1.7 x cm-^) and BiaSeg gOl (- 
10^^ cm“^). The defects are labeled by comparison with the 
previous results El within the close-up picture of Fig.|^b). 
The two types of Sb vacancies (Vsbi, Vsb 2 ) exhibit depres¬ 
sions of different lateral size at positive sample voltage, while 
the two antisite defects (Sbxe, Tesb) appear as protrusions of 
different lateral size. Obviously, the acceptors (Vsbi, Vsb 2 , 
SbTc) outnumber the donors (Tesb) by far as found consis¬ 
tently on all probed areas of the sample. This is in line with 
the known p-type conductivity of Sb 2 Te 3 . 

To this end, we used the identification of different defects in 
order to determine defect densities for several 40 x 40 nm^ 
surface areas as summarized in table U For two of these ar¬ 
eas, we additionally performed FF spectroscopy in order to 
determine E-p, — Ep. The FF spectroscopy recorded at dif¬ 
ferent spatial positions and a single field (6.7 T) within the 40 
X 40 nm^ areas revealed that Ep, shifts by less than 5 meV, 
which points to more long range potential fluctuations. Note, 
however, that the spatial sensitivity of LLq for potential fluctu¬ 
ations is limited by the magnetic length Ib = \jh/ (eB) being 
Ib = 10 nm at B — 6.7 T 041 . This also prohibits to detect 
potential fluctuations caused by individual defects, since the 
LL wave functions always cover several defects. 

From the counted amounts of different defects, we determined 
straightforwardly the overall defect density np by dividing the 
sum of all defects by the probed area. A first estimate of the 
effective p-type doping results from subtracting the num¬ 
ber of n-type defects from the sum of the p-type defects and 
dividing by the product of probed area and thickness of a QL 
(1 nm). This assumes singly charged defects due to the fact 
that their charge is not known ll54l and cannot be determined 
by STS. The numbers for Ep, np, and pes are also displayed 
in table lijrevealing, e.g., average values of no = 0.035 nm“^ 
andpgff = 3.1x 10^® cm“^. The hole concentration of Sb 2 Te 3 
has previously been deduced from transport measurements to 
be 5.3 • 10^® cm“^ for thin films grown by MBE ll58l and to 
be 2.7 • 10^® cm“^ for a polycrystal |j59l, both fits reasonably 
with the effective p-type doping deduced from our STM im¬ 


ages. 

By comparing area 1 and area 2, one observes that an in¬ 
creased peff leads to a larger value of Ep — Ep as ex¬ 
pected. Moreover, the increase of Ep — Ep (^ 25 %) is 
smaller than the increase in p-type doping (^ 40 %) indi¬ 
cating an increasing density of states at Ep with increasing 
distance from Ep, as expected for the Dirac cone. Quantita¬ 
tively, the expected charge carrier density of the Dirac cone 
of n{E) = {E — Ep)^/ {inh^Vp) amounts to n(166 meV) = 
2.8 • 10^^ cm“^, respectively, n(204 meV) = 4.2 • 10^^ cm“^, 
which nicely fits to the defect density observed by STM (ta¬ 
ble]^ when assuming one single charge per defect. However, 
this excellent agreement is probably coincidental, since, on 
the one hand, STM is only sensitive to the first QL Il43l . while 
DFT calculations show that the hole part of the TSS penetrates 
by more than one QL into the bulk of the substrate 14911^ . 
which increases the effective doping with respect to the one 
deduced from the counting of defects in STM measurements. 
On the other hand, the lower part of the TSS energetically 
overlaps with bulk states also to be doped by defects, which 
effectively decreases the achievable Ep—Ep at a given dopant 
density. 

The fluctuation of the number of defects in different areas 
of 40 X 40 nm^ is remarkably large. For Vsb, we find an 
average number of A^Vsb = 45 with a standard deviation 
of CTVsb = 10.5, while the statistical fluctuation should be 
a = 6.7 only. The deviation between crvsb ^ I’Y nearly 
60 % is significant regarding the six sample areas probed. The 
same trend, but without statistical significance, is also found 
for the Sbxe defects. Thus, we believe that the potential fluctu¬ 
ations are not only given by a statistical distribution of defects, 
but are also influenced by the kinetics of defect formation de¬ 
pending critically on local temperature na and local flux of 
the different constituents. This implies stronger potential fluc¬ 
tuations on large length scales as anticipated by a statistical 
distribution of acceptors only ll63l . 


















5 



sample voltage (mV) 




Energy (meV) 


FIG. 4. (color online) (a) dI/dV{V) spectra recorded at B = 0 — 7T as marked on a different sample area as the LLs from Fig. [^a) 
(14tab = 0.3 V, /stab = 100 pA, Vmod ~ 2 mV). For direct comparison of the spectra, each spectrum (at a distinct field), which is again an 
average of 10 spectra, is measured on the same point away from the defects. The vertical dotted line indicates the field-independent zeroth 
Landau level at the Dirac point Eo ■ The spectra are shifted vertically for clarity. Inset: dl/dV{V) spectrum for a larger energy range measured 
at iJ = 0 T in the same sample area (14tab = 0.3 V, /stab = 100 pA, Vmod = 2 mV), (b) Landau level energies deduced from Lorentzian 
fits of the curves in (a) and plotted versus momentum as deduced from eq.|^ The dashed lines are linear fits to the lower and upper part, 
respectively, with resulting Dirac velocities dd marked, (c) Full width at half maximum for each Landau level, as deduced by a Lorentzian fit 
function, plotted as a function of energy. Different B fields are indicated. The dashed line is a linear fit forced through AE = 0 meV at Ep. 


V. DIFFEIfENT DIRAC VELOCITIES EOR HOLES AND 
ELECTRONS 

Figure j^a) shows the LL measurement for area 2 at vary¬ 
ing B. At B > 2 T, quantization peaks are visible, which 
become more pronounced with increasing B. The field- 
independent LLq, identified as the Dirac point Ed, is found 
at Ed = 204 meV (see table |^. Its field-independence corrob¬ 
orates again the Dirac fermion nature of the TSS. In area 2, 
Ed could not be identified at B = 0 T due to a low dl/dV 
intensity of the TSS within the band gap (inset in Fig. Qa)), 
which is most probably due to a different density of states of 
the microtip in comparison with the microtip used to probe 
area 1 (inset of Fig. [^a)). This underlines the advantage of 
LL spectroscopy to determine Ed precisely. 

Figurej^b) shows the derived Landau level energies En, again 
deduced from the peak values of Lorentzian functions fit¬ 
ting the dl/dV curves. They are displayed with respect to 
their momentum kn as derived from the Onsager semiclassi- 
cal quantization condition ll64l . This condition provides an 
equation for the fc-space area Sn of the nth Landau orbits, 
namely S'„ = 27 r(nH- 7 )(l/(^), with 7 = 0 being the phase off¬ 
set for Dirac fermions 1651 . Close to Ed, the Dirac cone in 
Sb 2 Te 3 is isotropic ia such that Sn = leads to: 

The corresponding En{kn) exhibits the known linear band 
dispersion, however, with different ud for the lower and the 
upper part of the Dirac cone, namely wq = 4.1 ± 0.1 • 10® m/s 
and z;d = 4.9 ± 0.2 • 10® m/s, respectively. A slower ud 


below Ed has also been found for the 3DTI Bi 2 Se 3 1^ . 
The value for the lower part of Sb 2 Te 3 (0001) reason¬ 
ably fits with ARPES measurements of the same crystal 
(ud = 3.8 ± 0.2 • 10® m/s) iHa . The results also roughly 
agree with the ud deduced from DFT calculations, which 
reveal vd = 3.3 • 10® m/s for the lower and vd = 5.0 • 10® m/s 
for the upper part of the Dirac cone within the same momen¬ 
tum range ll4^ . Note, that a ud of 4.4 ± 0.1 • 10® m/s is 
obtained if evaluated in the same way as in Fig. |^b), which 
nicely agrees with vd = 4.2 ± 0.1 • 10® m/s deduced in 
Fig.|^b) and ud = 4.3 • 10® m/s found in ||43]| . 


Notice, that the area encircled by a full cyclotron orbit, as 
semiclassically required for Landau quantization lIMIl . con¬ 
tains a significant amount of defects Nd- For example, we ob¬ 
serve the LL with n = 1 on the electron side at i? = 2 T (Fig. 
I^a)) indicating the presence of a fully developed cyclotron 
orbit. Taking the cyclotron radius Vc = y/2n + 1 - Is =32 
nm, we get Nd = it ■ ■ ud = 130 defects. DFT shows that 

the penetration depth of the TSS on the electron side is about 
1 QL 1491 . such that this number is indeed a good estimate 
for the number of defects encircled. In order to set this ap¬ 
parently large number into perspective, we estimate the cross 
section for electron-defect scattering classically assuming the 
full cyclotron orbit and get CT 2 d < ('\/2 • 27rrc • = 0.8 

A in 2D, or regarding the depth of the QL (1 nm) a reasonable 

3D cross section of a^D < 8 A . 
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VI. QUASIPARTICLE LIFETIMES 

Finally, the energy dependence of the quasiparticle lifetime 
is extracted from the peak widths AE of the LLs as deduced 
from the Lorentzian fits. The result is shown in Fig. [^c) 
revealing a nearly linear increase of AE with energy. Recall 
that we operate at an energy resolution 6E = 4 meV being 
smaller than the determined peak widths. A linear fit to the 
data in Fig. [Jb), forced to AE — 0 meV at Ep, results 
in AE = (0.084 ± 0.003) • {E - Ep) (dashed line). In 
addition, one could anticipate a slight decrease of AE around 
i^D. We will offer a different interpretation of the apparent 
jump at E' ~ 250 meV later, but would like to stress that the 
linear behavior below Ed remains, even when excluding the 
data above Ed. The dip and the linear slope was observed 
similarly by LL spectroscopy on Sb 2 Te 3 thin films ll43l and 
on Bi 2 Se 3 bulk material cleaved in-situ im. Interestingly, 
the linear slope of AE(E) towards Ep is nearly identical 
for the Sb 2 Te 3 thin film, if applied for E — Ep > 50 meV 
(~ 0.085) ll43l . and only about 30 % smaller for the Bi 2 Se 3 
bulk crystal (~ 0.06) HD, albeit, in both cases, the dip at Ed 
is more pronounced. This is even more remarkable, since the 
defect density of the Sb 2 Te 3 thin film is a factor of two lower 
than the defect density Ed of the bulk crystal studied here. 
The peak width being larger than SE could be firstly caused 
by disorder broadening. Basically, electronic states in B field 
have a width of about Tc and are mostly localized by the 
disorder ll^ WT\ . Consequently, states with different origin 
in space and, thus, different potential energy spatially overlap 
resulting in a finite width of the peak in dl / dV, which, in first 
order, depends linearly on the product of potential gradient 
and Tc ll68ll . Since Tc increases with Landau level number n, 
the disorder broadening must also increase with |n|, which 
is in contrast to the experimental result on the hole side of 
Ed. However, it might account for the dip of AE(E) at Ed. 
The effect is most likely less pronounced in our experiment, 
since we operate at a larger Ed — Ep than for the thin films 
BD, such that the linear contribution of AE(Ed — Ep), 
being of different origin, is already larger. In line, within 
areas of about 40 x 40 nm^, {r^ = 40 nm, e.g. at E = 3.5 T 
and n = 4), the peak position of LLq does not shift by more 
than 5 meV, which is on the order of the scattering of AE, 
showing that disorder broadening is not the dominant effect 
for AE. 

Thus, we attribute the peak width to inelastic scattering, 
which could be either driven by electron-electron or by 
electron-phonon interaction. The energy dependence of 
electron-phonon interaction is typically dominated by the 
available phase space for electron scattering 16^ |70l . at 
least, at energies above the optical phonon energy being 
Elo = 20 meV for Sb 2 Te 3 E] E^. For LL distances 
smaller than Epo^ one could anticipate the Dirac cone density 
of states E(E, E = 0 T) c>c |E — Ed| as the available phase 
space, while for larger LL separations the scattering should 
be restricted to one Landau level exhibiting a degeneracy 
Dia,{E,B) oc E. Thus, one would expect AE to increase 
with |n|, i.e., away from Ed for |n| > 1 and for the well 
separated LLq, one would expect AE cx E. At energies 


below Ed, bulk states start to overlap with the Dirac cone 
II 49 I - I 52 I increasing the phase space for electron scattering 
even further, such that one would expect shorter lifetimes on 
the hole side than on the electron side in obvious discrepancy 
to the experiment. Also the other trends anticipated for 
electron-phonon scattering are not observed experimentally. 
Consequently, by the exclusion principle, we believe that 
electron-electron scattering, which typically implies an 
increasing scattering rate with distance from Ep GH is 
dominating AE. Indeed, at E = 0 T, one expects, a linear 
energy dependence for a Dirac cone with Ed = Ep Gl or 
for a doped Dirac cone at E > Ed — Ep being interpolated 
by a more quadratic dependence towards Ep GSl. The 
numbers calculated for graphene and graphite in these studies 
GIGS are surprisingly similar to the ones we observe for 
Sb2Te3(0001), e.g., being AE = 20 meV at E — Ep = 200 
meV Gl- However, this might be a coincidence, since 
a quantitative comparison would require a more detailed 
calculation including the E field, the penetration depth of 
the surface state into the bulk, which effectively changes the 
dielectric constant, and the bulk valence bands as additional 
scattering channels. Notice that the peak widths jump from 
about 15meV to about 22meV at an energy of ~ 250 meV. 
This coincides with the onset of an increasing dl / dV signal 
(Fig. I^a)) most likely being related to the onset of the bulk 
conduction band. This onset is also found in two-photon 
ARPES data to be about lOOmeV above the Dirac point 1761. 
Thus, interband scattering sets in as soon as the Landau levels 
from the Dirac cone overlap with the bulk conduction band. 
The fact that time-resolved ARPES data at E = 0 T find re¬ 
laxation times of about t = 1 ps at Ed — Ep = 0 — 100 meV 
G61G3, which would correspond to AE = h/r = 0.7 meV 
indicates that either the E field and the related localization 
of electrons or the presence of the STM tip might change the 
scattering rates significantly. 

Notice that the lifetimes observed by the fit of the line 
widths of ARPES data i?!] GSl are typically even shorter 
than the ones observed by LL spectroscopy MM and, 
thus, in stronger discrepancy to the time-resolved data 
lIMl ISTIl . This is most likely due to disorder averaging, i.e. 
bands from areas with different Ed (see table [I]i overlap 
in the E(fc) representation of the ARPES data. Notice 
further that a dominating influence of electron-electron 
interaction as a relaxation channel of hot electrons has 
also been deduced from the decay of standing waves from 
step edges of Bi 2 Se 3 ( 0001 ) thin films probed by STM at 
E = 0 T ll82l . The authors deduced that r decreases with 
|E —Ep| finding r = 100—10 fs at |E —Ep| = 0.05 — 0.3 eV. 


VII. SUMMARY 

In summary, we probed the Landau quantization of the 
topological surface state of Sb 2 Te 3 ( 0001 ) by STS for varying 
magnetic fields perpendicular to the sample surface. Our data 
reveal the Dirac type energy dependence, E„ cx -^InjE, with 
different Dirac velocities for the hole and the electron branch. 
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The Dirac point energy deduced from the zeroth Landau level 
spatially fluctuates by about 40 meV which could be traced 
back to varying densities of p-type defects beyond a statistical 

distribution. The peak width of the Landau levels is found to VIII. ACKNOWLEDGMENT 

be /S.E ~ 0.08 • \E — E-p \ pointing to a dominating influence 

of electron-electron interaction. We gratefully acknowledge provision of the sample by M. 

Wuttig, technical help by M. Pratzer, and financial support by 
the DFG via SFB 917, project A3 and Mo 858/13-1 as well as 
of Fonds National de la Recherche (Luxembourg). 
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